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INTRODUCTION 


Water droplets and ice crystals, the constituents of clouds, are very 
nearly transparent (i.e. they absorb almost no radiation in the visual wave- 
lengths (400 nm to 700 nm) ) . Clouds are also optically thick with optical 
depths for a one kilometer path ranging from 10 to 50, depending upon droplet 
size and number density of droplets. Therefore, visible light which enters 
a cloud is scattered many times before being absorbed or exiting the cloud. 
This type of process is well described by a diffusion model. 

In this paper it is shown how the radiative transfer equation reduces 
to the diffusion equation. To keep the mathematics as simple as possible, 
the approximation is applied to a cylindrical cloud of radius R and height h. 
The diffusion equation separates in cylindrical coordinates and, in a sample 
calculation, the solution is evaluated for a range of cloud radii with qloud 
heights of 0.5 km and 1.0 km. 

The simplicity of the method and the speed with which solutions are 
obtained give it potential as a tool with which to study the effects of 
finite-sized clouds on the albedo of the earth-atmosphere system. 


THE DIFFUSION APPROXIMATION 


The diffusion approximation has long been used in nuclear reactor theory 
(refs. 1-3) and has recently been applied to the transfer of visual radiation 
in snow (ref. 4) . In the diffusion approximation, the radiation is assumed 
to have an almost isotropic angular distribution, so that the specific 
intensity at space point £ for radiation traveling in direction ^ is 


I(r, fi) J(r) - 3Dfi • ^J(r^) 
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is the mean intensity, D is the diffusion coefficient, and V. is the gradient 
operator. Because the net vector flux 


$ = / M dfi 


(3) 


is proportional to the gradient of J : 
$ = - 4ir D V J 


it follows that J satisfies the diffusion equation 

2 -2 
V J - L J = 0 


(4) 


(5) 


where the diffusion length, L, is related to D according to 
= D/k 

K is the absorption coefficient. 

The net flux through a surface with normal h is 

$ = h . $ = _ 4 t7p 

n - 8n 


(6) 


where 3/3n denotes the directional derivative. The flux in the direction of 
+h is 


$ = ttJ - 2TTD 3J/3n (8a) 

n 


and the flux in direction -h is 


$ = ttJ + 2irD 3J/3n (8b) 

n 


The plus and minus signs indicate the positive and negative senses of h. 
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The diffusion coefficient and the diffusion length are related to the 
extinction coefficient x and the phase function of the cloud droplets 
according to 


D = l/[3x(l - (Jjg)] 


(9) 


and 


= l/[3x(l - o>) (1 - cog)] 


( 10 ) 


In these equations, co is the single scattering albedo (or fraction of light 
scattered in a single interaction with a cloud droplet) and g is the mean 
cosine of single scattering. Equations (9) and (10) are essentially approxi- 
mations to the dispersion relation derived by Mika (ref. 5) for the largest 
singular eigenvalue in the singular eigenfunction solution for a plane- 
parallel atmosphere and expanded by van de Hulst (refs. 6 and 7) in his 
scaling laws. It may be noted that the diffusion equation (5) is identical 
to the Eddington approximation in a plane-parallel medium, in that for this 
geometry, the second moment of the radiation field, K, in a conservative 
atmosphere is equal to one-third of the mean intensity 

K 5 j ;^1 I dy == j J (11) 

deep in the medicmi. Equation (11) uses the assxmption of the plane-parallel 
atmosphere and describes the direction of propagation ^ in terms of a polar 
coordinate system with a polar angle 6 = cos“^y such that 0 = 0 is perpendi- 
cular to the plane of symmetry. 


THE CYLINDRICAL CLOUD 


Clouds in the atmosphere vary markedly in shape and size. To keep the 
mathematics as simple as possible, we have applied the diffusion approximation 
to a cylindrical cloud of radius R and height h. The cloud is illuminated 
from the top by a diffuse source, which is normalized to unit flux. It is 
assumed that none of the radiation that escapes the cloud returns. 

With this choice of geometry and boundary conditions, the diffusion 
equation (5) separates. The solution for the mean intensity may be found 
as an expansion of standard mathematical functions in the form 


J(r, z) 


CO 

Z [a exp (a z) + B exp (-a z)]j„(3 r) 
-m m m m Om 

m=0 


( 12 ) 


1087 



with JQ(gjj^r) being the zeroth order Bessel function of the first kind. The 
separation constants and gui related through the diffusion length in 
the following manner; 


2 £>2 ^ ^-2 

a = g + L 
m m 


(13) 


3m is the solution of the transcendenta l equation that results from the appli- 
cation of the boundary condition at the sides; 


e R J (3 R) - C J„(3 R) = 0 
m 1 m 0 m 


(14) 


where 



(15) 


Equation (14) is solved numerically using a Newton-Raphson method of root 
finding. The application of the boundary conditions at the top and bottom 
yields expressions for the expansion coefficients 


^ l+2Da 

j _ 2 m 

m it3RJ,(3R) 

m 1 m [(l+2Da ) -(l-2Da ) exp (-2a h) J 

mm m 


[l+( 



(16) 


and 


A 

m 


(2Da -1) 
m 

(2Da +1) 
m 


exp (-2a h) 
m 


B 

m 


(17) 


Our primary interest in this study is in the fate of the energy incident 
on the cloud. As one might expect, the energy may be reflected back out the 
top, may be transmitted out the bottom, escape out the sides, or be absorbed 
within the cloud. The power incident on the cloud top is 


T-.'t' /•2ir ,, ,R .(“) , 2 

^ (z=0)dr = 1 TR 

TOP 0 0 z 


( 18 ) 
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If we normalize the reflected and absorbed power by this amount, we find that 
they may be expressed in terms of the expansion coefficients and geometrical 
properties of the cloud in the form 

E ' = 2ir E [A (l+2Da ) + B (l-2Da )]J, (g R)/8 R (19) 

TOP „m m m mlm m 

m=0 


= 2ir E [A (exp (a h)-l) + B (1-exp (-a h))] 
SIDE „ m ^ m m ^ m 

m=0 


[J (B R) + 2DB J (B R)]/a R (20) 

0 m mlmm 


e'*' = 2t7 E [a exp(a h) (l-2Da ) + B exp (-a h) 

BOTTOM „ m m mm m 

m=0 

(l+2Da )] J, (B R)/B R (21) 

mlmm 

and 


^ (exp(a h)-l) + B (1-exp (-a h) )] • 

ABS ^ m m m m 

m=0 


J (B R)/a B R 
1 m mm 


A SAMPLE CALCULATION 


( 22 ) 


In a sample calculation the above expressions were evaluated for clouds 
of heights 0.5 and 1.0 km and a range of cloud radii from 0.5 to 10.0 km. 

The cloud droplet radius was assumed to be 10 microns with a number density 
of 100 cm”^. The mean cosine g was assumed to be 0.8516 and the single 
scattering albedo ai was chosen to represent nearly conservative scattering 
(1-0) = 10”®) and non-conservative scattering (l-o) = 10 ”^). These values are 
typical of cumulus clouds as seen in Deirmendjian's model Cl (ref. 8). 

The effect of the finite radius is marked, as can be seen in Figs. 1 and 

2. The albedo of an isolated cloud, 1 km thick, may be reduced by 5 percent 

or more if the radius is less than 5 km. This reduction in albedo is due to 

I 

the leakage of energy out the sides of the cloud. As the radius becomes ' 
smaller, the escape out the sides becomes more and more important. Clouds 
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whose radii are about equal to their height lose nearly as much energy out the 
sides as is reflected back out the top. 

As the radius increases we expect the results to closely approach those 
from a plane-parallel treatment. Table 1 shows the difference between the 
results using the diffusion approximation for a cylindrical cloud and those 
using the Eddington approximation for a plane-parallel layer. The difference 
is quite small for clouds whose horizontal extent is much larger than the 
vertical . 

It is generally known that Monte Carlo techniques, which are currently 
being used to study the effects of finite-sized clouds (refs. 9 and 10), 
consume great amounts of computer time. The above results using the diffusion 
approximation required less than 30 seconds execution time on a CDC CYBER 175 
computer. The simplicity of the approximation and the speed with which 
results are obtained give the diffusion approximation potential as a tool 
to study the effects of finite-sized clouds on the earth-atmosphere system. 


CONCLUDING REMARKS 


Clouds represent an optically thick medixam for visible radiation in which 
the internal radiation field is very nearly isotropic. Such a medium is well- 
suited to a description by a diffusion model. Applying the diffusion 
approximation to a cloud of cylindrical geometry, the fraction of the incident 
energy emerging from each of the cloud's surfaces has been calculated. The 
amount of radiation escaping from the sides becomes significant when the 
cloud's horizontal extent is less than ten times its vertical extent. The 
speed and simplicity of the method argue for its use to study the effects 
of finite-sized clouds on the earth's albedo. 
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TABLE I . -COMPARISON BETWEEN EDDINGTON (PLANE-PARALLEL) AND 
DIFFUSION (CYLINDRICAL) APPROXIMATIONS 
(Cloud Height = 1 km) 



Plane 

-Parallel 



Cylindrical 






R = 10 km 

5 km 

2 km 

1 km 

0.5 km 

-ft 

Albedo 

0.875 

0.848 

0.821 

0.744 

0.631 

0.464 

e 

u 

O 









Transmission 

0.125 

0.114 

0.104 

0.076 

0.040 

0.009 

-2 

Albedo 

0.549 

0.538 

0.528 

0.498 

0.450 

0.361 

1-W =10 








• 

Transmission 

0.009 

0.009 

0.008 

0.007 

0.005 

0.002 
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Figure 1.- Fate of energy incident upon a cylindrical cloud of 
height 0.5 km as a function of radius (solid line: l*-w « 10~^; 

broken line: 1-w = 10~2) , 


FRACTION 

OF 

INCIDENT 

RADIATION 



Figure 2.- Fate of energy incident upon a cylindrical cloud of 

height 1.0 km as a function of radius (solid line: l-o) = 10“^; 

broken line: l-w = 10”2). 



